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Abstract 

A method to extract .B-physics parameters (6-quark mass and fs, 
Jb s decay constants) from currently available lattice data is presented 
and tested. The approach is based on the idea of constructing ap- 
propriate ratios of heavy-light meson masses and decay constants, re- 
spectively, possessing a precisely known static limit, and evaluating 
them at various pairs of heavy quark masses around the charm. Via 
a smooth interpolation in the heavy quark mass from the easily ac- 
cessible charm region to the asymptotic point, -B-physics parameters 
are computed with a few percent (statistical + systematic) error using 
recently produced Nf = 2 maximally twisted Wilson fermions data. 

1 Introduction 

Heavy flavour physics is a corner of the Standard Model where chances are 
higher to uncover signals of new physics [TJ [2]. However, to extract from 
experiments useful phenomenological information, it is mandatory to have an 
accurate knowledge of the relevant hadronic matrix elements of the effective 
weak Hamiltonian. For low mass states (up to around the charm mass) lattice 
QCD (LQCD) represents the ideal framework where such calculations can be 
performed with well under control systematic errors [3]. 

Due to present day computer limitations, it is not possible, however, 
to work directly with the heaviest quarks (as the 6-quark) propagating on 
the simulated lattice. Various strategies, more or less inspired to the heavy 
quark effective theory (HQET) [U [5] , have been devised to circumvent this 
intrinsic difficulty, which go from non-perturbative matching of HQET onto 
QCD [5] to finite size scaling methods with relativistic heavy quark(s) [TJ. 
Relativistic heavy-quark actions designed (highly tuned) to have reduced 
cutoff effects [8] have also been employed for this purpose. Encouraging 
results have been obtained by several groups [9] though different is the level 
at which the various relevant systematic effects are controlled. 

In this work we wish to present a novel approach to -B-physics in which 
the fe-mass point is attained by interpolating from the charm region to the 
asymptotic infinite mass regime suitable ratios of heavy-light (hi) meson 
masses and decay constants, computed at a number of pairs of quark mass 
values lying slightly below and somewhat above the charm mass. The key 
feature of the approach is the use of ratios of physical quantities which by 
construction have a well defined and exactly known infinite /z-quark mass 
limit. Injecting knowledge of meson masses and decay constants at the charm 
region, their /i-quark mass evolution can be computed by a chain of successive 
steps up to values as large as about twice the charm mass. The 6-physics 
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region is finally reached through an interpolation from the simulated points 
to the exactly known infinite /i-quark mass value. 

A first test of the viability of the method is presented here. It has been 
carried out by exploiting the unquenched Nf = 2 data recently produced 
by the ETM Collaboration [TOj E] which makes use of maximally twisted 
Wilson fermions [T2] . The results obtained in this feasibility study are very 
encouraging and compare nicely with the unquenched determinations today 
available in the literature [13] as well as with PDG numbers pp. We get 



where, as indicated explicitly in eq. ( II. ip . the b- mass has been run in a world 
with two (active) flavours. The results in eqs. (II. ip to (11.31) represent a 
"first principle" determinations of i?-physics parameters with errors whose 
magnitude can be systematically reduced. The quoted uncertainty will be 
discussed in sects. 2 and 3 for eqs. (11.11) and (ll.2p - (ll.3l) . respectively. 

A few observations are in order here. First of all we would like to remark 
that the results above are extracted from unquenched LQCD data where 
u and d light fermions are dynamical, while heavier quarks are introduced 
only as valence quarks. This scheme is what goes under the name of "par- 
tially quenched" setting (see ref. [H] for a discussion within the twisted mass 
regularization of QCD). Systematic errors due to partial quenching are not 
included in the figures quoted in eqs. (11.11) to (II. 3p . The second observation 
is that no complicated renormalization steps are required for the method to 
work, because, as noted above, the necessary inputs are (ratios of) physical 
quantities (M-pseudoscalar meson masses or decay constants) evaluated at 
/i-quark masses around the charm region which are extracted from the ex- 
isting (large volume) lattice configurations produced for the study of pion 
physics. 

The central value of fi^ IS ' Nf ~ 2 (fi^ IS ' Nf ~ 2 ) in eq. (II. ip may look somewhat 
higher (though still compatible within statistical errors) than the available 
phenomenological estimates of the MS 6-quark mass at its own scale, which 
lie in the range 4.2-4.3 GeV. However, it is not unlikely that, when the 
quenching of quarks heavier than u and d will be removed with the inclusion 
of dynamical s and, possibly, c quarks, nib will receive corrections which 

one can argue will tend to make the quantity /if f,S ' 7V/_4 (/i^' S ' Ar/_4 ) somewhat 
smaller than the number given in eq. (11.11) p. 

1 Indeed, if we evolve the intermediate result p,^ IS ' Nf ~ 2 (2 GeV) = 5.35(32) GeV from 



„ MS,N f =2 , . MS,N f =2 
^b \^b 



) = 4.63(27) GeV, 



(1.1) 
(1.2) 
(1.3) 



f B = 194(16) MeV, 
f Bs = 235(12) MeV, 
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The results (11.21) and (11.31) . instead, are only affected at a level of less 
than 1% by our present uncertainties in the 6-quark mass, because the hu/d- 
and /is-meson decay constants happen to have a rather mild dependence on 
the /i-quark mass. 

The content of this paper is as follows. In sect. |2]we discuss the theoretical 
basis underlying the strategy that we propose to extract the value of the b- 
quark mass from present day LQCD data and we provide a rather accurate 
determination of it with controlled errors. In sect. [3] we extend the method to 
the determination of the Jb and fs a decay constants. We conclude in sect. |4] 
with a few words on how to improve the quality of the numbers (II. ip to ( II. 3p . 
and how to extend the present method to other /i£-physics quantities the large 
/i-quark mass behaviour of which is known. We defer to an Appendix some 
technical details concerning the way chiral and continuum extrapolations of 
hi meson masses and decay constants are performed. 



2 6-quark mass 

In this section we present a simple strategy aimed at determining the value 
of the 6-quark mass through a smooth interpolation of suitable ratios of hi 
pseudoscalar lattice meson masses from the well accessible charm region to 
the asymptotic (infinite mass) point where these quantities have an exactly 
known value. Inspired by HQET results we consider the lattice ratios (a = 
lattice spacing) 



y L {x^\X-fi^a) = 



n = 2,---,N. (2.1) 



In eq. ( 12. ip and in the following by a "hat" we denote quark masses renormal- 
ized at 2 GeV in the MS scheme. By fi£ we indicate the renormalized light 
quark mass, while > {x^ ^ are pairs of (renormalized) "heavy" valence 
quark masses lying around (from below to somewhat above) the charm mass. 
The function p(log fih) is the factor that "transforms" the renormalized MS 
quark mass at 2 GeV scale into the so-called "quark pole mass" . In formulae 

p(log£ft)Afc = /i£ ole - ( 2 - 2 ) 

In continuum perturbation theory (PT) p is known up to N 3 LL (i.e. up to 
next-to-next-to- next-leading-log) order terms included [THl US] . Finally TV is 



2 GeV to the 6-quark mass scale by using anomalous dimension and /3-function of the 
Nf — 4 (rather than Nt = 2) theory, the value of the 6-quark mass gets lowered by about 
3% compared to the value we give in (jl.ll) . 
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the number of /i-quark masses at which the values of the hi pseudoscalar 
lattice meson masses, Mfc e , are supposed to have been measured. 

The choice of the form of eq. (12. ip is suggested by the HQET (continuum) 
asymptotic equation [U [5] 



lim — -M- = constant 7^ . (2.3) 

pole . ..P°^ c 



Although the above constant is known to be 1, its value is not really needed 
here. 

In order to simplify our subsequent analysis we keep fixed the ratio be- 
tween two successive values of the heavy quark masses in eq. (12.11) . Calling 
it A > 1, we set 



A % {n-l) ,,(""!) x (n) ' ~(n) ' 



4 n) ^ ^ 

Notice that in p^/p^ ^ the mass renormalization constant factor, Zp 1 , 
cancels out q 

Ratios of the kind defined in (12.11) are introduced with the idea that they 
might have a smoother chiral (fie — > ^ u /di with fi u / d the light quark mass 
that yields the physical value of the pion mass) and continuum limit than 
each of the individual factors. Setting 

y(x {n \ A; pu/d) = lim lim y L {x {n \ A; p e , a) = 

! M hu/d (l/xW) p(logAxW) 
M hu/d (l/\xW) pQagxW) ' 1 ' ' 

where we have introduced the (continuum limit) shorthand notation 

M hu / d (l/x) = M hu/d (l/x,p u / d ) , 

we observe that (for all A > 1) eqs. (12. 2p and (12. 3p imply the following exact 
property 

\imy(x,X;p u/d ) = 1. (2.6) 

From lattice data the function y(x, A; p u /d) can be determined at certain 
discrete values of x (x^ n \ n = 2, • • • , N). In order to extend our knowledge 



2 We recall that in maximally twisted LQCD the twisted mass renormalizes according 
to // = Zp l {i. If standard Wilson fermions were to be employed, the quantity rhh = 
Zg (rriQh — m CI ), should be used in place of in eq. (|2.1[) . 
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outside these particular points, while at the same time fully exploiting the 
strong constraining power provided by eq. (12.61) . we imagine proceeding in the 
following way. Suppose the perturbative expansion of p has been computed 
and resummed up to N P LL order. Then we can define a tower of y-ratios, 
y\v> V — 0, 1, • • • , P + 1, such that 

y(x,X;p u/d )\ -1 0( 1 ), (2.7) 
P v (logx) p+i/ 

provided p in eq. (12. ip is correspondingly taken at tree-level in the case of 
p = 0, or to N P_1 LL order for p > 0. Then for sufficiently small values of x 
we parameterize y\ p in the form 

y(x, A; pL u/d ) = 1 + 771 (log x, X; p u / d )x + i] 2 (hgx, A; p u / d )x 2 , (2.8) 

where the coefficients r)j(logx, A; p u /d), 3 — 1 3 2, are p-dependent (though 
to lighten the notation we do not display explicitly this dependence in the 
following), smooth functions of logx which tend to zero as A — > 1 and to 
some fixed constant [U [5] . With the ansatz f)2.8p and at any order 

where PT results for p are available, it is not difficult to determine the i]j 
coefficients from lattice data, assuming that their log x-dependence can be 
ignored in the range of masses where the above formulae are used. 

It is important to remark that the ansatz (12.81) is based on the same kind 
of assumptions under which HQET is usually employed in the study of heavy 
quark physics. A posteriori, we check that the best fit values taken by the 
coefficient functions rjj come out of a reasonable order of magnitude. Indeed 
we find that i]ir and 772^ are 0(1) quantities. 



2.1 Implementing the method 

Let us start considering for concreteness the case where p is taken up to LL 
order and subsequently compare the results we get in this way with what one 
would obtain taking p at NLL-order or at tree- level (at tree- level p = 1). 

In order to determine the coefficient functions rjj we proceed as follows. 
Let us make for the smallest p,h value the choice p,^ = 1-230 GeV (we recall, 
we are referring to the MS scheme at the scale of 2 GeV). Fixing A = 1.278 
(see below) and, in this exploratory study, iV = 4, we shall successively 
consider the /i-quark masses 

= 1.230 GeV, 
/4 2) = A/4 1} = 1.572 GeV, 
pf = A 2 /^ = 2.009 GeV, 

= A 3 /^ = 2.568 GeV . (2.9) 
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a~\L 


, 3 x T) 


a/J,£ — a^, se a 




a/j, h 


3.80 


24 3 


x 48 


0.0060, 0.0080 
0.0110, 0.0165 


0.0200, 0.0250 
0.0300, 0.0360 


0.2700, 0.3100 
0.3550, 0.4350 
0.5200 


3.90 


24 3 


x 48 


0.0040, 0.0064 
0.0085, 0.0100 
0.0150 


0.0220, 0.0270 
0.0320 


0.2500, 0.3200 
0.3900, 0.4600 


3.90 


32 3 


x 64 


0.0030, 0.0040 


0.0220, 0.0270 


0.2500, 0.3200 


4.05 


32 3 


x 64 


0.0030, 0.0060 
0.0080, 0.0120 


0.0150, 0.0180 
0.0220, 0.0260 


0.2000, 0.2300 
0.2600, 0.3150 



Table 1: Lattice size, light (= sea), strange- and charm-like bare quark mass 
values used in the analysis presented in this work. The number of correlator 
measurements was 240 in all cases, but for (3 = 4.05, where it was 130. The 
r /a-values 4.46(3), 5.22(2) and 6.61(3) are employed at f3 = 3.8, (3 = 3.9 
and (3 = 4.05, respectively [T 7 ^ . The overall scale and the light quark mass are 
set by the experimental values of f n and m n via chiral fits of the pseudoscalar 
meson mass and decay constant data in the light quark sector FTS\ \T^jj. Here 
we use fi u /d = 3.6(3) MeV and r = 0.433(14) fm. 

Actually at each lattice spacing we will be dealing with the dimensionless 
quantities /tjfVo, j = 1,2,3,4, of which the numbers quoted in eq. (12. 9p 
represent the central values in physical units. Uncertainties on r$/a and Zp l 
(present at the level of about 3%) will be taken into account in the final error 
analysis. 

From the set of the ETMC simulation data [TU] with parameters detailed 
in Table (TJ we extract the values of the hi pseudoscalar meson masses that 
correspond to the fih values listed in ( I2.9p . With these masses we construct 
the lattice ratios (12.11) on which a combined continuum and chiral fit is per- 
formed. As we hoped, ratios appear to have a mild dependence on the light 
quark mass fi£ and small cutoff effects, as seen for instance in fig. [TJ This 
makes our continuum and chiral fit straightforward and numerically robust. 

The red squares in fig. [2] represent the numbers = y{x <y1l \ 1.278; p> u /d)\i, 
n = 2,3,4, computed at the = 1//4™'' values in the list (12. 9p with the 
ratios p(log/t| l ?l ~ 1 ^)/p(log/i| l ri ^) computed at LL order (i.e. p — 1 in eq. (12.71) ). 
The best fit through the red squares and the point at x = determines the 
values of the rjj coefficients (j = 1, 2) and yields the middle (red) curve in 
the figure. We note that a second order polynomial in x is necessary to get 
a good fit to the data (a straight line forced to pass through the point y = 1 
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Figure 1: 



Lattice 



spacing 



Mfo(M 3) ; fa, a)/M^(fa^; fa, al- 



and fit dependence of the ratio 
The black square with its error is the 
combined continuum and chirally (fie — > fi u /d) extrapolated value. Here and 
in all the following figures uncertainties possibly affecting the value of the 
variable in the horizontal axis are propagated to the quantity plotted on the 
vertical axis. 



at x = would have a very large x 2 )- The quadratic fit gives for the quan- 
tities r r]i and r§f]i numbers of order unity, in agreement with the standard 
assumptions underlying HQET. 

At this stage, having in our hands the quantities y[ n ^ 
for any n (actually for any x), the iterative formula 



y(xW, 1.278; fi u/d )\, 



(2) (3) 



•Vi 



(K+l) 



A 



-K 



M hu /d(P>h 



M hu / d (fr h 



0-h 



p(logAk 



P^Ogflh , 



;2.io) 



should be looked at as a relation between the mass of the hu/d- meson, 
Mhu/d(f'h <:+1 ^), and the corresponding heavy quark mass fi\^ +1 \ which is 
fully explicit if the initial, triggering value Mh U /d{n) ls assigned. The lat- 
ter can be accurately measured, as fij^ lies in the well accessible charm quark 
mass region. We show in fig. [3] the quality of the continuum and chiral ex- 
trapolation of the triggering mass lattice data. Once this number is known, 
determination of the 6-quark mass is tantamount to find the value of K at 
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Figure 2: Continuum data for y\o (blue dots), y\i (red squares), y\2 (green 
triangles). The corresponding best fit curves are drawn with A = 1.273 (lower 
curve, in blue) and A = 1.278 (middle curve, in red and upper curve, in 
green). In all cases fi£ —> //«/<£. 



which Mhu/difrh ) takes the experimental B- meson mass value, Mb- Call- 
ing Kb the solution of the resulting eq. (12.101) (as shown in fig. HJ we find 
Kf, = 6), one gets for fi^ the simple formula (valid for renormalized as well 
as bare masses) 

h = ^ Kh ti ] - (2-11) 

A few related remarks are important here. 1) It is not really necessary 
to have the lattice M pseudoscalar meson masses computed at values of fit 
matched exactly as indicated in eq. (12.91) . A \ih interpolation between nearby 
Mfc e masses can be carried out if necessary. This is what we have actually 
done in the numerical study we present in this paper. 2) It is not a priori 
guaranteed that eq. (12.101) can be solved for an integer value of the exponent 
K. This is not a problem, however, as one can always retune the parameter A 
(and at the same time readjust the values in the sequence (12.91) ). so as to end 
up with an integer for Kf, (this is the reason why the peculiar value A = 1.278 
was chosen). Alternatively one could adjust the starting value of the heavy 
quark mass or both. 3) A detailed discussion of the numerical analysis will be 
given in a forthcoming publication [19]. Here we only mention that a simple 
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jjf, ~ 1.23 GeV 
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3.5 
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Figure 3: The M^{(x h \ fie, a) lattice data extracted from the simulations de- 
tailed in TableUi The black square with its error is the continuum and chirally 
extrapolated value giving M hu / d {fi~^) = 1.89(10) GeV. 



SU(Nf = 2) chiral perturbation theory NLO-formula was used to model the 
fie dependence of the triggering hi meson mass and y-ratios, while 0(a 2 ) 
effects have been parameterized (at each fi h ) by fie independent terms. A 
few further details on this point are given in Appendix A. 

Following the procedure outlined above, one finds the result given in 
eq. (11. II) . or equivalently the renormalization group invariant (RGI) value 

/if GI,Nf=2 = 7.6(5) GeV, (2.12) 

where in the running only two flavours are assumed to be active and the 
conventions of ref. [20] for the RGI quark mass have been used. 

2.2 Discussion and error budget 

It is important to check the degree of reliability of the key smoothness as- 
sumption we have been implicitly making on the function y(x) and test the 
sensitivity of the procedure and its result (12. lip to the order of PT at which 
the expansion of p is truncated. To this end we have repeated the entire anal- 
ysis above using for p both a lower (tree-level) and a higher (NLL) order per- 
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Figure 4: The numerical solution of eq. A2.10\) giving K b = 6. 



turbative approximation in place of the previously employed LL order trunca- 
tion. We recall that in the large /i-quark mass limit y\o = y(x, A; fi u /d)\o and 
2/ 1 2 = y(x, X; fi u /d)\2 approach 1 (see eq. (12. 7p ) with corrections 0(l/logx) 
and 0(l/(loga;) 3 ), respectively. 

One finds a (very) little shift in the value of A necessary for the solution 
of the corresponding iterative equation ( I2.10p to be an integer (from 1.278 to 
1.273) if we go from y\i to y|o- The shift is instead totally negligible within 
our statistical accuracy if we move from y\i to y\2. 

The fits to y\o and y\2 data are shown in fig. [2] (lower blue and upper 
green curve) together with the fit to y\x data (middle red curve). One clearly 
sees that the y\ p curves tend to become flatter and flatter as we move from 
p = to p = 2. As for the values of the 6-quark mass, instead, the number 
extracted from 7/|o-ratios is only about 2% smaller than the one obtained 
using the y|i-ratios. The difference between the latter and the one extracted 
from £/| 2 -ratios is smaller than 1%. 

The stability of the value of the 6-quark mass with varying p should not 
come as a surprise. It is enough to notice that, if the 7/| p -ratios were exactly 
known, for a generic value of p one would get (recall eq. (I2.10p ) 



[K+ih 



p(logAk 



■yiH 3) ---y ( P K+1) 
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,(2)„,(3) „,(A'+1) _ x -K M hu/d{P h 



yf'yf" ■ ■ ■ y { K+1) = *~ K > ^ 



as all the intermediate p factors (except the first and the last) cancel out in 
the l.h.s. leaving behind simply the product of the y|o-ratios. The small p 
dependence we have found in the value of A (hence in //&) is due to the slightly 
different level of accuracy by which the 7/| p -ratios (which instead significantly 
depend on p) can be described by a polynomial in x with the lowest order 
coefficient set to unity. In this respect, increasing p is expected to improve 
the quality of the ansatz (12.81) and reduce the systematic error associated to 
it. This is so until p becomes so "large" that the accuracy of the p estimate 
gets spoiled by the renormalon ambiguity in its perturbative expansion [21] . 

To account for the truncation to LL of the perturbative expansion for 
p we have conservatively decided to attribute to the 6-quark mass value a 
systematic error of 1%, which is added in quadrature to the other errors 
discussed below, leading to the total error quoted in eq. (11 .ip (and (I2.12p ). 



2.2.1 Error budget 

The total error we attribute to the 6-quark mass results (11. ip and (I2.12p takes 
into account a number of statistical and systematic effects which we now 
briefly illustrate. The relative error on the product of the continuum ^/-ratios 
in the l.h.s. of eq. (I2.1(jp is only about 1%, whereas the pseudoscalar meson 
mass in the charm region (Mwd^ ) in the r.h.s. of eq. (I2.10p ) contributes 
a relative error of about 5%. These errors are the result of our statistically 
limited knowledge of hi- meson correlators, r^/a and Zp d as well as of a 
number of further systematic errors. Among the latter we mention those 
coming from the fit ansatz underlying the combined continuum and chiral 
(pi — > p u /d) extrapolation, the error due to the x-interpolation to the b- 
mass point, as well as the (tiny) error inherent the numerical solution of 
eq. (12.101) (giving Kb = 6 and A = 1.278). As we discussed above, the 
effect on pb due to the truncation of the p perturbative series to order p 
is very small, not larger than 1%. Another .5-1% systematic error comes 
from the possible (neglected) logarithmic dependence of the r]j, j = 1,2 
coefficients. The relative uncertainty on rjj associated with these effects can 
be estimated to be 0(a s (l/x)) ~ 10 — 15%, a number which is never larger 
than the statistical errors on their best fit values. Finally cross-correlations 
between the different quantities (stemming from common ensembles of gauge 

3 At the moment the statistical error on Zp quoted by ETMC is about 3% at the 
simulated lattice spacings. 
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configurations) are as usual taken into account by a bootstrap error analysis. 
Further technical aspects of the error analysis are deferred to ref . [TH] . 

The information provided in figs. [TJ [2] and [3] about the a 2 , fie and x 
dependence of the intermediate quantities entering our analysis as well as 
about the precision in solving eq. (I2.10p (see fig. HJ) shows that the global 
systematic uncertainty is well within (or below) our present statistical errors. 

We conclude by observing that, as expected, our results for the 6-quark 
mass (and fs or /# s discussed in the next section) do not significantly depend 
on the value of the intermediate quantity tq which is only employed to ease 
continuum extrapolations, while the physical scale is ultimately set by /„.. 

2.3 b and c quark masses 

Although not necessary, the phenomenological value of the D-meson mass 
could have been used as a triggering mass. In this would have to 

be identified with \i c . 

We note in this context that, since we get (see the black square in fig. [3] 
Mhu/difah) = 1-89(10) GeV, i.e. a number that practically coincides with 
the experimental value of Mo (M D o = 1.865 GeV pQ), our method imme- 
diately yields for the charm mass the estimate fi c = 1.23(06) GeV. For the 
phenomenologically important b- over c-mass ratio we then get 

^ = 4.31(24), (2.14) 

in very good agreement with other estimates @. 

In closing this section we note that an independent determination of the 
6-quark mass can be obtained repeating the same analysis as before but us- 
ing Mb s , instead of Mb, and replacing fi u / d with fi s . By doing that we find 
a result which is fully consistent with the one in eq. fl2.12p . Alternatively, 
and perhaps more interestingly, one could use Kf, as determined from Mb to 
predict Mb s , or better the ratio MbJMb, by the method we are proposing 
in this paper. Such an analysis is in progress and will be presented else- 
where |19j . 

4 It is interesting to compare our unquenched result for the ratio of RGI masses 
I,Ni ~ 2 /m? I ' Nf ~ 2 = 0.232(13) (the inverse of eq. (I2.14jl ) with the correspond- 
ing quantity [m? GI ' Nf=0 = 1.654(45) GcV]/[mf GI - N f= Q = 6.758(86) GeV] = 0.245(7) 
determined using the quenched data of ref. [55] (for the c-mass) and [23] (for the 
6-mass). More recently the work of ref. |24] has appeared where the number 

Tficf ' Nf ~ 4 (3 GeV)/m^ IS ' Ni ~ 4 (10 GeV) = 0.273(3) is quoted. This value is well consis- 
tent with our result (12TT41) which translates into m c ' f (3 GeV)/m* ■ f (10 GeV) = 
0.274(15). 
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3 fs and fs s decay constants 



A strategy very similar to the one outlined in sect. [2] can be employed to 
extract accurate values of the fs and /b s decay constants from available 
lattice data. In analogy with what we have done before, one should now take 

, \ , x m/2 fhijl/x) CfQogjxX)) [p(logx)] 1 / 2 

1 ' '™ Ul/xX)' Cr(logx) [ P (logAx)]V2 V- 1 * 

with the (continuum limit) shorthand notation 

fhi(l/x) = f M (l/x,fi e ) , 

where the quark mass ji£ must be extrapolated to either fi u / d or to the ap- 
propriate strange quark mass value, jl s , depending on whether one wants to 
compute f B or f Bs . 



1.10 




Figure 5: Lattice spacing and p,£ dependence of the ratio 
f^ e (jl^-, fif ,a)/ f^ e (jl^] fie,a). The black square with its error is the 
combined continuum and chirally (p,g — > £i u /d) extrapolated value. 

The form of the function z(x, A; jig) is dictated by the continuum asymp- 
totic formula 



lim 

x->0 



'p(logx) fheO-/x) 
x C^ at (\ogx) 



constant ^ 



(3.2) 
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which follows by matching HQET to QCD [U [5]. The presence of the factor 
(jstat comes f rom the fact that in HQET the axial (and vector) current needs 
to be renormalized. The renormalization constant C^ at is known in PT up 
to three loops [25]. The ratio of p factors (raised to the appropriate power) 
is there to convert MS heavy quark masses to pole masses (see eq. (12.21) ). 

As before, the function (13. ip has been defined so as to fulfill the exact 
asymptotic constraint 

lim z(x, A; fx A = 1 , (3.3) 

x— >0 

from which the small x expansion 

z(x,X,fx t ) = 1 + Ci(logx, X;fxe)x + ( 2 (logx, X;fi e )x 2 , (3.4) 

follows. Again the coefficients Q (log x, X; fie), j = 1,2, are smooth functions 
of log x which tend to zero as A — > 1 . 

In analogy with what we did in sect. |2]in determining the 6-quark mass, 
with the purpose of checking the robustness of the procedure, we shall take 
(jstat an( j p a £ i ncreas i n g orders in PT, from tree-level up to NLL order, and 
construct z^-ratios endowed with the asymptotic behaviour 

z(x,X;fi u/d ) -1 0(- — l —— -) , (3.5) 
p v (logx) p+1/ 

Just like in the case of the determination of the 6-quark mass, the values of fg 
or fs s that we shall extract will be almost independent of the PT truncation 
order. 



3.1 Implementing the method. The case of fs 



The ^-ratios (13. ip have been evaluated at the reference /i-quark masses of the 
list (12.91) for each of the lattice spacings and light quark mass values given 
in Table ([T]). When we perform the continuum and chiral extrapolation of 
the ETMC lattice data for the ratios (13.11) of hi pseudoscalar meson decay 
constants (again based on simple chiral NLO-formulae supplemented with fif 
independent 0(a 2 ) corrections - see Appendix A), as hoped, a rather smooth 
behaviour is found since most of the a 2 and fit dependence gets canceled in 
taking the ratio. The observable dependence on fi£ and a 2 is mild and/or 
hardly significant within our present statistical errors (see e.g. fig. ED- 

From the structure of eq. (13 .ip one derives the iterative formulae (analo- 
gous to eq. flUOD with >) = z(x^,X; fi u /d)\ P ) 



7 (2).(3) . . . JK+1) 
^p 



(3.6) 



X K/2 



fht(P> 



h j 



sistat 



(l0| 



cr(iogfcf +i) )\\ p(io g/ ii ] 



p(lO[ 



ah 



p = 0,l,2 
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Figure 6: Continuum data for z\q (blue dots), z\\ (red squares), z\i (green 
triangles). The corresponding best fit curves are drawn with A = 1.273 (upper 
curve, in blue) and A = 1.278 (middle curve, in red and lower curve, in 
green). In all cases fi£ — > fi u /d- The blue and red vertical lines represent the 
location of the b-mass as extracted from y\ andy\\ data (with A = 1.273 and 
A = 1.278,), respectively. The green vertical line is practically on top of the 
red line and it is not visible. 

Similarly to what we did in fig. [21 we collect in fig. [6] continuum and chirally 
extrapolated data for z\ p , p = 0, 1, 2 and best fit curves through these data 
and the value at x = 0. Thus, for instance, the middle (red) curve is the 
parabola (eq. fl3.4[) ) which best fits the values of z^ = z(x^ n \ 1.278; jj, u /d)\i, 
n = 2,3,4, at the heavy quark masses (12. 9p . The red vertical line marks 
the position Xf, which corresponds to the previously determined value of /2& 
(eq. (I2.12p ) and crosses the curve at the point z[ Kb ^ = 1.278; fr u /d)\i- 
With the help of this number and the values of z± for 4 < j < + 1, 
eq. (13. 6p provides a determination of fhu/difib) in terms of fhu/difih ) (with 
LL-accurate fit for the z-ratios). As observed before, the latter does not 
necessarily has to be identified with the phenomenological value of fo for 
the method to work, as what we actually need to know is the dependence of 
fhu/difah) on frh at around the charm mass. Nevertheless, since, as remarked 
in sect. 12.21 M/ lu /rf(/i[ l 1 ' ) ) coincides with the experimental value of Md, we are 
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in position of evaluating f D , obtaining f D = fhu/di.^c) = 211(9) MeV, com- 
patible with the result fr> = 197(9) MeV given in ref. [IT]. The latter was 
obtained in the standard way (see ref. [10]) from the same ETMC gauge con- 
figuration ensembles, but with a rather different analysis method where the 
meson masses rather than the renormalized quark masses were kept fixed as 
a — > 0, resulting in somewhat different statistical (no use of Z P ) and system- 
atic errors as compared to the present study. For the present computation of 
fhu/difth ^) = Id the quality of the continuum and chiral extrapolation of our 
lattice data is shown in fig. [71 Taking as triggering value the continuum and 
chirally extrapolated value of the pseudoscalar decay constant computed at 

-(1) A- 

f B = hu/d{h) = 194(16)MeV , (3.7) 
which is precisely the result (11.21) . 




Figure 7: The f^ e (fi h , jlf ,a) lattice data extracted from the simulations de- 
tailed in TableUi The black square with its error is the continuum and chirally 
extrapolated value giving fh u /d(f l h\ Au/d) — 211(9) MeV. 



3.2 Discussion and Jb error budget 

To test the reliability of the interpolation of our trial functions, z(x), to 
6-mass point, we have explicitly checked the stability of the result (13.71) to 
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the order of PT at which the expansion of p 1//2 /C^ ta * is truncated. For this 
purpose we have repeated the whole previous analysis employing values of 
p 1 / 2 '/C^ at computed at tree-level (p = 0) and NLL (p = 2) order. Upon 
comparing with the decay constant values obtained in these other ways we see 
that numbers obtained using the z|o-ratios (upper blue curve in fig. [6]) differ 
by less than 1% from the value one gets from the red data (LL z |i-ratios). 
If one employs the lower green data (coming from the NLL z| 2 -ratios) the 
difference with the previous determination is totally negligible (about 0.1%). 
This specific systematic effect on fg was hence conservatively estimated to 
be ~ 0.5% and added quadratically to the full error. 

As in the case of the determination of the 6-quark mass, the remarkable 
numerical stability of fs with varying p can be traced back to the good 
quality of the interpolation ansatz (13.41) and the relation (again valid for 
exactly known 2| p -ratios) 



Q^(iogP4 J 



. M) (3) ...JK+1) 

^p ^p ^p 



(2) (3) JK+1) _ xK /2 fhe(fj>h i > 



3.2.1 The fs error budget 

The total error we attribute to fs in eqs fll.2p and (13. 7p comes in almost equal 
parts from the product of ^-ratios in the l.h.s. of eq. (13. 6p and the value of 
fhu/d(fah ) an d is a combination of statistical and systematic errors stemming 
from the same sources already illustrated in the case of the 6-quark mass in 
sect. (12.21) . As we saw above, the systematic error stemming from the trunca- 
tion of the PT expansion of p 1 / 2 /C^ at has a negligible impact on Another 
1-2% systematic uncertainty comes from the possible (neglected) logarithmic 
dependence of the Q, j = 1, 2 coefficients. The relative uncertainty on Q as- 
sociated with these effects can be estimated to be 0(as(l/x)) ~ 10 — 15%, 
a number which, as in the case of the i]/s entering our analysis for is 
never larger than the statistical errors on their best fit values. In any case 
inspection of figs. [5] to [7| shows that all systematic uncertainties are smaller 
than our current statistical errors. 



3.3 The case of f Bs 

In order to come up with a determination of fs s one has simply to re- 
peat the whole procedure setting p^ — > p s . With reference to the value 
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of fif s ' Nf - 2 (2 GeV) = 99(7) MeV given in [55], one finds (thanks to the 
equality of pL^ with the charm quark mass) 

f Ds = hsiffl) = 252(7) MeV (3.9) 

and the best fit z-ratio curves shown in fig. [HJ We may quote as our final 
result (see also eq. (11. 2j ) 

Sb. = hs{M = 235(12) MeV , (3.10) 

where errors are estimated as in the case of . We also note that the result 
(13.91) for the decay constant of the D s meson is in agreement with that of 
ref. [llj and contributes to further reduce the possible "tension" with the 
recent Cleo data reanalysis [27]. 

A more complete analysis of the many possibilities and refinements one 
can envisage will be presented in |TS]. 

1.18 
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Figure 8: Same as fig. m — > fj, s . 

As is clear by comparing the results for f B and fg 3 (eqs. f !3.Tj) and (I3.10p . 
respectively), our method yields a significantly smaller error for the decay 
constant of the strange S-meson. The reason is that no large statistical fluc- 
tuations from the light (u/d) quark propagators nor (valence) chiral extrapo- 
lation uncertainties enter the computation of fs s ■ In view of this observation 
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we remark that, if one would know with high accuracy the ratio /b//b s , a 
more precise determination of fs could be obtained by simply multiplying 
this number by fs a El Actually the quantity fs / fB s can be accurately com- 
puted by a simple generalization of the method discussed in this paper. It is 
indeed sufficient to consider the double ratio 

w{x,X;fi u/d ,fi s ) = [fh/ud/ fhs]{^/x)[f hs / f hu/d }{l/x\) = (3.11) 

and follow the procedure we described before starting from the triggering 
quantity [fhu/d/fhs\iPh)- This kind of analysis is under way and will be 
discussed elsewhere [T9] . 

4 Conclusions and outlook 

In this paper we have proposed a novel strategy to determine 5-physics 
parameters from currently available Wilson fermion simulation data. As 
a first test of the method, we have computed in (the continuum limit of) 
QCD with Nf = 2 light dynamical quarks the (renormalized) 6-quark mass 
as well as the 5-meson decay constants, fs and fs s , employing the gauge 
configurations recently produced by the ETM Collaboration with maximally 
twisted Wilson fermion action [10J. 

The method provides rather accurate numbers with errors that are domi- 
nated by the uncertainties related to the limited statistical accuracy by which 
the (two-point) hi pseudoscalar meson correlators and the quark mass renor- 
malization constant, Zp 1 , are evaluated. A better assessment of the system- 
atic errors due to the limited knowledge of logarithmic corrections can only 
come from data taken at quark masses larger than the ones displayed in 
eq. O. 

In several respects the present feasibility study could benefit from the 
nice properties (particularly 0(a) improvement [29KH]) of maximally twisted 
Wilson fermions. Indeed, an important feature of the present computation 
is the pretty good control we have of cutoff effects, which (judging from the 
spread between values at the coarsest lattice spacing and those at the con- 
tinuum limit) are always smaller than 10%. This is so both for the triggering 
quantities at the charm mass scale and for the y- and ^-ratios, which involve 
higher quark masses (up to twice the charm mass). 

It is also interesting to note that the whole procedure only relies on the use 
of physical quantities that can be easily determined from lattice simulations, 
while the need for a renormalization step is limited to establishing the relation 

5 This observation is far from original, see e.g. ref. [28,. 
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between the renormalized charm-like mass and the values of the triggering 
pseudoscalar meson mass. Fixing this relation requires the knowledge of Z P . 
No extra renormalization factor is needed for the calculation of the decay 
constants of interest if maximally twisted fermions are used as the charged 
axial currents are exactly conserved at finite lattice spacing. 

There is a lot of room for improvement in the application of the method, 
like reducing the statistical error of the correlation functions, using several, 
suitably smeared meson sources, increasing the accuracy by which Zp is 
known and incorporating in the analysis the new ETMC set of data that are 
coming out at a finer lattice spacing (j3 = 4.2). 

Needless to say, the method can be straightforwardly extended to LQCD 
computations with u, d, s and possibly c dynamical quarks where quenching 
uncertainties are virtually absent. In this respect we wish to note that in 
simulations with Nf = 3 dynamical quarks, although low energy hadronic 
effects in the S-meson wave function are correctly treated, a conflict remains 
between the number (Nf = 3) of dynamical quarks running in the loops 
(and thus relevant for the subtraction of UV divergencies) and the number 
(Nf = 4) that instead should be used for continuum RG-evolution at scales 
above, say, 1.5 GeV. This problem and the related RG-uncertainties are 
completely removed if also the c quark is made dynamical. 

Finally, we remark that the strategy we have outlined can be applied to 
any other hi physical quantity the large \ih behaviour of which is known 
(typically from large quark mass arguments). 

Acknowledgments - We wish to thank R. Sommer and N. Tantalo for 
useful discussions and all the other members of ETMC for their interest in 
this work and a most enjoyable and fruitful collaboration. 



Appendix A — Chiral and continuum extrapolations 

For the reader's convenience we collect in this Appendix the standard formu- 
lae we have used to perform the necessary chiral (fii — > (i u /d) and continuum 
(a — > 0) extrapolations of our lattice data on pseudoscalar masses and decay 
constants. 
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Fit ansatz for hi meson masses and decay constants 

We have modeled the /^-dependence of Mm and fht in a form that is con- 
sistent with the known results of (NLO) SU(2) chiral effective theories for 
pseudoscalar mesons made up by a light plus a heavy quark. The form of 
the fit ansatz was chosen generic enough so as to encompass (see discussion 
below) the expected /^-dependence both in the case when the heavy quark 
is treated as static J3UJ ED E2] and when the latter is considered non-light 
but still relativistic [331 El] ■ 

In ro-units we write for masses and decay constants 

a 2 

M M r = Co + C 1 ft e r + —C L , (A.l) 

f M r = D + D^ro + d x ^% log + %D L , (A.2) 

where the last terms, with parameters Cl and Dl, have been included to 
cope with the expected 0(a 2 ) discretization effects. The fit parameters C...'s 
and D 's in general depend on /t^, though the form of the fit ansatz flA.ip 
and (IA.2P has been actually employed only for a fixed value of the heavy 
quark mass, namely for fih = fi^ \ when we evaluate the so-called triggering 
meson mass and decay constant, respectively. 

We checked that, within the statistical accuracy of our data, no fii depen- 
dence is visible in the 0(a 2 ) terms in eqs. flA.ip and (IA.2I) . The fit ansatz for 



Mh£r does not include logarithmic terms. This is consistent with the results 
of the chiral effective theory for hi pseudoscalar mesons with a light plus a 
non-light and relativistic quark [331 El], but it can be equally well regarded 
as a simple Taylor expansion leading to a polynomial interpolation of data 
points with a very smooth dependence on £i e (see fig. ED- In this sense the fit 
ansatz for M^r is also consistent with the spirit of the effective theory for 
static-light mesons [301EDE2], where the fi£ dependence of M M is expected 
to be a tiny effect (as we indeed find). 

The coefficient d\, multiplying the term ~ /^log(/t^) in eq. flA.2j) . was 



taken as a free fit parameter. Numerically we get for di/D® at fih = fi^ a 
value (—1.0 ±0.4) which falls in between (and agrees within statistical errors 
with) what is expected from the arguments of refs. [301 ED E2] (where the 
heavy quark is treated as a static source) and those of ref. [33] (where it is 
taken as a relativistic particle). The result we find for d\ is not surprising as 
our heavy quark mass lies in the charm region. 

The low energy constant fo and Bq have been taken from recent ETMC 
analyses of light meson quantities [T7j . 
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Fit ansatz for ratios 



For the quantities 

^(/4r } ;/^ a ) /m(a! w) ;A^q) n = 2 3 4 ( A 3 ) 

M w(/4 n ~ 1} ; A/» a) ' /w(Al n_1) ; a) ' 

which enter the ratios y and 2 at the various /^-values, we employed fit 
ansatz analogous to (and derived from) eqs. (IA.10 and (IA.20 above. They all 
are of the following form: a leading term plus a term linear in £i£ and another 
one proportional to a 2 (and fa- independent) . 

Note that, if we assume that d\ does not appreciably vary as fa changes 
by the factor A ~ 1.27-1.28 (which is a natural expectation in any effective 
theory for hi mesons), one finds that in the ratios (IA.3I) the possible fa log(fa) 
dependence cancels at NLO, and is pushed to NNLO. 
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